Dynamics of the Formation of Bright Sohtary Waves of Bose-Einstein Condensates in 

Optical Lattices. 
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We present a detailed description of the formation of bright sohtary waves in optical lattices. To 
this end, we have considered a ring lattice geometry with large radius. In this case, the ring shape 
does not have a relevant effect in the local dynamics of the condensate, while offering a realistic set up 
to implement experiments with conditions usually not available with linear lattices (in particular, 
to study collisions). Our numerical results suggest that the condensate radiation is the relevant 
dissipative process in the relaxation towards a self-trapped solution. We show that the source of 
dissipation can be attributed to the presence of higher order dispersion terms in the effective mass 
approach. In addition, we demonstrate that the stability of the solitary solutions is linked with 
particular values of the width of the wavepacket in the reciprocal space. Our study suggests that 
these critical widths for stability depend on the geometry of the energy band, but are independent 
of the condensate parameters (momentum, atom number, etc.). Finally, the non-solitonic nature of 
the solitary waves is evidenced showing their instability under collisions. 



I. INTRODUCTION 

One of the most relevant characteristics of Bose- 
Einstein condensates (BEC) is the coherent nature of 
their macroscopic wavepackets 0. This is particularly 
important in the case of BEC in optical lattices, where 
the continuous diffraction of the wave in the lattice dis- 
torts severely the dynamics expected for a free conden- 
sate. 

The behavior of BEC in optical lattices has drawn the 
attention of the research community almost since the first 
realizations of condensates in the laboratories . Several 
experiments have demonstrated fundamental aspects of 
the lattice dynamics as Bloch oscillations |1| as well as 
the progressive refinement of the techniques of loading 
and control of the condensate inside the lattice (J]. 

The affinity between the non-linear Schrodinger equa- 
tion (NSLE), which governs the mean field dynamics of 
the condensate, and the equation for the propagation of a 
light pulse in the slowly varying envelope approximation, 
suggests that several phenomena of non-linear optics can 
be reproduced in BEC. Among them, the formation of 
self-trapped (i.e. non dispersive) pulses is especially 
relevant, due to its applications in signal transmission. A 
fundamental property of the NLSE is its integrability 6] 
and, as a consequence, the existence of soliton solutions. 
Bright solitons, or more generally solitary waves, in BEC 
have been demonstrated experimentally for the case of 
attractive interactions |3Ja and dark solitons in the case 
of repulsive interactions |9lll0j|. 

Bright solitons can be also obtained for repulsive inter- 
actions in condensates near the edge of the Brillouin zone 
of an optical lattice (gap solitons) 0, I IS] . The first 
experimental observation of this type of solitary waves 
has been reported very recently l3| . Despite the interest 
of this matter, little theoretical effort has been addressed 
to understand the mechanism of spontaneous formation 
of the solitary wave from an initial wavepacket. Only 



very recently, some studies include matter radiation as 
a fundamental process accompanying solitary waves |l3j | . 
It is the aim of this paper to clarify the essential points of 
this process with the aid of ab initio numerical solutions 
of the NLSE equation. Our calculations suggest that 
matter radiation is the fundamental relaxation process 
that transforms the initial wave to a solitary wavepacket. 
Accordingly to |0| , the origin of the radiation is the dis- 
parity of the wavefunction with the exact solution of the 
Schrodinger equation. In this paper, we further clarify 
this point by identifying a close link between the radi- 
ation process and the presence of high-order dispersive 
terms in the effective wave equation for the pulse enve- 
lope. These terms arise from the non-parabolic character 
of the energy band in the neighborhood of the packet's 
mean momentum. The effect of higher order non-linear 
terms has been studied in the context of light propagating 
through optical fibers as a source of non-integrability 
of the NSLE. 

The numerical results will be interpreted with the help 
of an analysis based in the Wannier function picture. 
Wannier functions have been employed before connecting 
BEC d yna mics in optical lattices and the Bose-Hubbard 
model [l6j as well as for obtaining approximated models 
(see for instance ^|) . However, it has been rarely used 
as a tool to understand the dynamical aspects of BEC in 
optical lattices as provided by the exact numerical inte- 
gration of the NLSE. Our approach will be summarized 
in section ^ where the effective equation for the enve- 
lope of the wavepacket in the Wannier basis is derived, 
including the higher order dispersive terms usually ne- 
glected in soliton theory. We shall use this approach to 
understand the results of the integration of the NLSE in 
section IlIII 

Our numerical computations are two-dimensional. For 
convenience, we shall use ring optical lattices instead of 
the usual linear ones. In addition to represent advantages 
for the numerical calculations, we believe that this geom- 
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etry is also practical in experiments where the structure 
of the lattice changes with position, or when different mo- 
menta have to be imprinted to wave packets at different 
positions. Up to the authors knowledge, these ring opti- 
cal lattices have not been yet created in the laboratories, 
although we find no fundamental difficulties in their con- 
struction (for instance in a similar way as proposed in this 
paper). We would like to emphasize that, since the ra- 
dius of the ring is taken much larger than the wavepacket 
dimensions, the dynamical effects related to the shape of 
the waveguide are small, therefore it is reasonable to ex- 
pect that the results described in this paper will also be 
applied in the case of other smooth geometries, as lin- 
ear lattices (however the discussion of collisions would 
require periodic boundary conditions). 

II. THEORY 

We consider a localized BEC released in an atomic 
waveguide at an initial time t = 0ms. Atomic waveg- 
uides of arbitrary shape have been created in laborato- 
ries by means of dipole-force confinement 0, |23| or 
current carrying structures . In the case of con- 

densates with low atomic density and strong transversal 
confinement, the dynamics of the wavepacket is essen- 
tially longitudinal, i.e. one-dimensional ■ In these 
circumstances, the transversal dynamics is frozen to the 
lowest energy eigenstate. Naming X2, a^s}, the co- 
ordinates along the transversal and axial directions of 
the waveguide respectively, the wavepacket factorizes as 
*(r,i) = £,i2{xi,X2)^3{x3,t). Where ^12(3^1, 2:2) corre- 
sponds to the lowest energy eigenstate of the transversal 
Hamiltonian. Therefore, the axial dynamics can be de- 
scribed by a one-dimensional Gross-Pitaevskii equation, 
which results from the projection of the original 3D on 
the transverse state, 

in^^Ux3,t) = {Ho+Hnl)U^3,t) (1) 

where Ho = -^SV^i ^(2:3), H^i = giD Mx3,t)\\ 

with giD — {Anh'^Nas/m) J J dxidx2\^i2{xi,X2)\'^, N 
being the number of atoms and the scattering length 
of the condensed element (in our case ^^i?6). The po- 
tential V{x3) corresponds to the axial lattice. Equation 
is a one-dimensional equation along the coordinate 
X3 , and ^3 (ccs , t) is normalized to unity. For simplicity, 
therefore, in the following we shall drop the coordinate 
index. 

In our study, the condensate is assumed to be released 
initially in a bare waveguide {V{x) = 0, at t = 0ms), 
with a particular momentum hko. Afterwards, the con- 
densate expands freely during 0.5ms, before the lattice 
potential is switched on linearly in time. After 10ms, 
the optical lattice reaches its maximum value, and re- 
mains constant, V{x) = Vbsin^(7rx/ao), for t > 10.5ms 
(ao being the intersite distance). Under these circum- 
stances, we have found that the condensate is loaded 



into the lattice adiabatically enough, with most of the 
population (typically about 85%) in the first band. In 
addition, the remaining population in the excited bands 
disperses rapidly over all momenta, reducing the ampli- 
tude probabilities accordingly. As a result, the relevant 
dynamics corresponds to the first band with little inter- 
action with the excited population. In the following, we 
shall assume ^{x, t) to represent only the wavefunction in 
the lowest band, normalized accordingly with its corre- 
sponding fraction of the total population. 

To understand the evolution of the condensate once 
the optical lattice has been generated, we find convenient 
to express the wavefunction either as a superposition of 
Wannier functions of the first band 

e(x,i)=^e*'V(^,Oa(x-^) (2) 
e 

with J2i 1/(^1 OP — 1; or the alternative superposition of 
Bloch functions 

ax,t)^ ydfce-"('=+'=«)*/''Cfc(i)V'fc+fco(^) (3) 

with Hoipk+koix) = e{k + ko)->jjk+koix)- These viewpoints 
are linked by the identity 

^k+M = -^^e'('=+'='')V (:.-£) (4) 

Using ||2Jl and ©, together with the orthogonality con- 
dition J dxa*{x — £)a(x — £') — 61^11, we find 

/(£,t) dfce"""=+'=°)*/''e*"Cfc(t) (5) 

Introducing Eqs. |(2J) and Q into and neglecting 
the overlap between Wannier functions of adjacent sites 
(assumption of tight binding) we find 

e^'-«'(^»?i^/(^,i))a(x-^)- 

y dfce-"('=+'=«)*/\(fc + ko)Ckmk+ko{x) + 

91D Y.e^'°' \a{x ~ a{x ~ e)\f{£,t)\^ f{e,t) (6) 
i 

In the cases discussed in this paper, the condensate 
wavepacket extends over some tens of lattice sites. Ac- 
cordingly, the corresponding wavefunction in the recip- 
rocal space is a narrow peak centered on fco, so we can 
approximate 

eik + ko) ~ e{ko) + hvgk + ^-^r - ask^ (7) 

where Vg = {l/h)(d/dk)e\k„ and 1/m* = 
{l/h'^){d'^ /d'^k)e\ko correspond to the standard def- 
initions of the group velocity and effective mass 
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respectively. Additionally we keep the higher order 
dispersion term as = — (l/3!)(9'^/9^A:)e|fco, which would 
be neglected if dealing with condensates extended over 
a greater number of lattice sites. 

Substituting Eq. and using and JSJ, the first 
term in the rhs of © can be written as 

-*«3^}mt) (8) 

Using this form and the orthogonality of Wannier func- 
tions, we can project Eq. © on the Wannier state 
e''°^a{x — tj to give 

|e(fco) + v.Hh)- - — — - ^a,—^^mt) 
+ ~g\fii,t)f fie,t) (9) 

where g = guj J dx\a{x — £)\'^. Defining adimensional 
space and time variables as 77 = {i—Vgt)/£o and r — t/ro, 
with io ~ y''hTo/\m*\ and tq a characteristic time, Eq. 
© can be rewritten as 

|q(?7,r)|%(7],r) (10) 

with 9(7;, r) = ^/|V^/(^7,r), a3\m*\^/yh'/X/\ 
and Eq = e(A:o)7o/?i. This equation contains a first cor- 
rection to the dispersion resulting from the non-parabolic 
shape of the band structure around fco- Since contains a 
third derivative it is expected to become important the 
wider the wavepacket is in the reciprocal space or, alter- 
natively, the closer is fco to the zero-dispersion point. In 
fact, the propagation of light in monomode fibers follows 
an equation, which is formally the complex conjugate of 
(Cnil JM- 

III. RESULTS AND DISCUSSION 

To test the above viewpoint and, particularly, the rel- 
evance of the higher dispersion terms in Eq. H10|) . we 
have computed ab initio numerical solutions of the Gross- 
Pitaevskii equation over a range of parameters. Sys- 
tematic calculations of the dynamics of propagating soli- 
tary waves in more than one dimension can be effectively 
speeded up with the choice of a waveguide geometry with 
periodic boundary conditions. On the other hand, our 
interest comprises also the study of solitary wave colli- 
sions. This kind of experiments require the formation of 



two solitary waves of different momenta, and represents 
a non-trivial extension of the present experimental se- 
tups used in the formation of bright gap solitons, which 
consider the translation of the whole lattice instead of 
imprinting an initial momentum to the wavepacket. To 
fulfill these requirements, without departing from a real- 
istic situation, we find advantages in proposing the use 
of ring-shaped optical lattices (see Fig. QJ. These kinds 
of lattices are composed by a ring-shaped atomic waveg- 
uide of the same sort as those used in [isL IT^ I20ll2lll2^ . 
The optical lattice itself may be generated by a periodic 
angular modulation in the transverse distribution of in- 
tensity of a laser (of incidence perpendicular to the lattice 
plane) using the appropriate mask. As in [23| the waveg- 
uide can be loaded in different locations using a properly 
shaped initial trap. Finally, different initial momenta can 
be communicated to wavepackets at every position by, for 
instance, a phase imprinting method jlOl l25l | . 

The ring-shaped atomic waveguide is defined by a po- 
tential 

K,g(r, (j), z) = ^moi^z^ + ]^mu}l{r - r^f (11) 

where in our calculations = 27r x 700Hz, ujr = 
2-K X 50Hz and ro = AO/jm. Since uj^ » (^r, for low- 
density condensates we can consider the dynamics along 
the z direction as frozen, and write in the same fash- 
ion as in section^ ^(r,t) = S,z{x)£^{r,(f),t), ^z{x) being 
the ground state of the harmonic oscillator in z direc- 
tion. Following the steps in section^ we find the two 
dimensional counterpart of Eq. 

^, , , /I d 1 32 V , , , 

(12) 

where g2D = {'iirfi^Nas/m) J dz\£^ziz)\*, and V{z) — 
Vb sin^(7rro(?!i/ao) is the lattice potential. We take as ini- 
tial time the moment in which the transference of the 
condensate to the waveguide has been accomplished. At 
t — Otos, an angular momentum L = h58 is imprinted 
to the condensate. This corresponds imprinting a linear 
momenta hko = L/tq along the waveguide circle. Note 
that, although high momentum rotating condensates are 
unstable in three dimensions and decay into smaller vor- 
tices, the waveguide reduces the effective dimensional- 
ity of the problem to one dimension. At i = 0.5ms 
the lattice potential is grown linearly in time during 
IOtos. In all the cases shown in this paper, the poten- 
tial has 128 wells and, therefore, the lattice constant is 
ao = 27rro/128 = 1.96/zm. From t = lO.Sms on, the 
lattice is stationary and the wavepacket evolves almost 
one dimensionally inside the waveguide, in the presence 
of the optical lattice 

All our numerical computations correspond to solu- 
tions of Eq. H12|) and, therefore, are two-dimensional. 



4 



However, since for our choice of ujr the radial size of the 
wavefunction in the waveguide is much smaller than the 
radius tq , the assumption of frozen radial dynamics made 
in section |n] is reasonable. This gives us the possibility 
of interpreting the results of the numerical integration of 
Eq. (|12|l with the theory developed in sectional with 
xi — z, X2 — r and X3 = r0 ~ ro<p- 

We have performed a set of numerical calculations in- 
volving a range of different lattice potential strengths Vq 
whose associated frequencies range from QAHz to 256Hz, 
number of atoms ranging from 50 to 200, and initial mo- 
menta hko ~ h58/rQ and hko — TiQA/rQ. To our knowl- 
edge, the experiments with condensates with these low 
number of atoms are difficult (for instance repeated mea- 
surements will destroy them ). However, the solitons ob- 
served recently fl^ are composed by 300 atoms. The 
scalability of the problem to higher number of atoms 
could be achieved provided the lattice intersite distance 
is preserved (therefore the band structure) and the initial 
condensate size is larger in X3 so that the atomic density 
at every site is kept similar (so that the tunneling rate 
is not changed). However, phase fluctuations may play 
an important role in these scaled systems and, therefore, 
could change drammatically the whole scenery and pre- 
vent the formation of the solitary wave. 

A typical process of the formation of bright BEC soli- 
tary wave is shown in Fig. |21 for a condensate of 100 
atoms. Figure |2a) shows the evolution of the envelope 
of the wavepacket projected into the Wannier functions 
of the first band {f{£,t) in Eq. ©). Since the initial 
momentum is close to the edge of the Brillouin zone, at 
t — Qms {Vq = 0) approximately 20% of the wavepacket 
extends into the second zone. As a result, even in the 
case of a perfect adiabatic switch on of the lattice, part 
of the population is loaded in the excited band. Our nu- 
merical results suggest that the population in the excited 
band expands rapidly in space, its probability amplitude 
decreasing accordingly, and therefore affects little to the 
evolution of the population in the lowest band, where the 
solitary wave is formed. The relevant dynamics, there- 
fore concerns the population of the first band. As an 
example of this, Figure[2Ib) shows the total atomic den- 
sity &i t — ISOtos, the solitary wave corresponds to 
the pulse of part (a), while the background noise corre- 
sponds to the radiated population in the first band and 
to the population of the excited bands. After the lattice 
is switched on, the condensate wavepacket in the first 
band is found to evolve freely with a velocity correspond- 
ing to the group velocity ~ Vg{ko). The stabilization of 
the pulse shape becomes apparent at times greater than 
t — 120tos. It is important to note the presence of a 
small trapped population near tq^ = 250/iTO in addition 
to the moving solitary wave. This population is emitted 
by the wavepacket during the time before adopting the 
solitary shape. Our calculations point out that this ra- 
diated population remains trapped in the lattice wells, 
since its low density prevents Josephson tunneling be- 
tween sites. This fact is supported by the inspection of 



Fig. El (a) , where the time evolution of the population 
in the first band, depicted already in Fig. |21 is plotted 
now projected on the Bloch basis. Also by the inspec- 
tion of Fig. O (b), where the same projection is plotted, 
but now restricted to the narrowest spatial box enclosing 
the solitary wave at every time. In both plots, the soli- 
tary wave corresponds to the peak near the edge of the 
Brillouin zone, and the smaller peak at zero momentum 
corresponds to the trapped radiation. In part (b), this 
later peak appears only at t = 120ms, i.e. when the ra- 
diation takes place, and disappears at later times, since 
the solitary wave moves out from the region where the 
radiation is trapped. The differences between the peaks 
belonging to the solitary wave in (a) and (b) suggest that 
there is an additional radiated component, which follows 
the solitary pulse (i.e. they have similar momenta), but 
located at a different spatial region. A close inspection of 
figure El demonstrates that it can be attributed to small 
amplitude fluctuations at the edge of the solitary wave. 
In the case of condensates of higher density, as in [Tsf . 
radiation trapping is not possible since its increased den- 
sity allows for Josephson tunneling. In these cases, the 
radiation peak is found at momenta different from 0. 

Figure Efa) and (b) also show that the solitary pulse 
is not located exactly at the edge of the Brillouin zone, 
but at some particular point inside. This is a remarkable 
point in the light of the discussion of the previous sec- 
tion, since implies the presence of non-vanishing third- 
order derivatives in the equation of motion l|10(l . Note 
also, that this also the case of any solitary wave with 
non-zero group velocity. It is known that the presence 
of a third derivative in the non-linear Schrodinger equa- 
tion prevents its integrability and, as a consequence, the 
appearance of solitons. In contrast, it has been reported 
the presence of solitary waves of light in optical fibers, 
which radiate constantly with a rate decreasing in time 
[Tsf . These sorts of solutions have a spectral distribution 
of the form shown in Fig. El 

From our calculations, therefore, it becomes apparent 
that radiation plays a fundamental role in the process 
of formation of solitary waves, since it represents a dis- 
sipative process that forces the convergence to a stable 
self-trapped solution. We find important to stress that, 
in contrast to other alternative schemes for bright soliton 
formation (for instance changing the sign of the non lin- 
ear parameter using Feshbach resonances), the presence 
of higher order dispersion is particular to the case of op- 
tical lattices, and provides a natural dissipative process 
leading to the convergence of the initial wavepacket into 
a solitary wave. 

The radiative process has also consequences in the dy- 
namics of the center of mass of the condensate wave 
packet, as a result of the conservation of the overall mo- 
mentum. This phenomenon is reflected in our calcula- 
tions, as it is shown in Fig. EI where the time evolu- 
tion of the momentum of the forming solitary wave is 
depicted as a function of time for two different initial 
momenta, hko = ?i58/ro and hko = fi-64/ro. It should be 
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stressed that the momenta are calculated from the mean 
value of the differential operator applied to the envelopes 
f[ro4',t) of the lowest band. Since a,i t — Oms part of 
the wavepacket extends beyond the limits of the first Bril- 
louin zone, the computed averaged momenta at this time 
is always smaller than fco- The first case (fco = 58/ro) 
corresponds to the same calculation as Fig. El the forma- 
tion of the solitary wave is described quite dramatically 
as a rapid stabilization of the wavepacket momentum. 
An examination of the time evolution of the wavepacket, 
shows that the presence of a strong oscillating region be- 
fore the solitary wave formation corresponds to the inter- 
action of the wavepacket with its radiation, as one runs 
over the other. In addition. Fig. 0] shows for this case 
an increase of the wavepacket momentum as the solitary 
wave is formed. As pointed out above, the low density 
of the radiated wave reduces drastically the Josephson 
tunneling rate. As a consequence, the momentum of the 
radiated wave is small, if not zero. In this situation, the 
momentum of the forming solitary wave should increase 
to account for momentum conservation, i.e. it will be 
accelerated. The case is rather different for the initial 
momentum hko = HQA/rQ, where the wavepacket mo- 
mentum decreases rather than increases. We have found 
this to occur in situations in which the initial wavepacket 
is closer to the edge of the Brillouin zone. In this circum- 
stance, the matter is radiated mostly to the excited band, 
and not to the lowest band. This effect can be attributed 
to Landau-Zcner tunneling induced by the nonlinear po- 
tential, which is more probable the smaller is the energy 
gap (i.e. at the edge of the Brillouin zone). 

As a general trend, we have found that the width of the 
condensate wavepacket in the reciprocal space is a critical 
parameter for stability. Figure El shows the time evolu- 
tion of this quantity for condensates with different initial 
number of atoms and for different initial momenta. In all 
these cases the formation of the solitary wave has its sig- 
nature in the stabilization of the width of the wavepacket 
to a rather similar value. The same behavior is also 
present for other lattice potentials, however the conver- 
gence resulting to a different value. This result seems 
to be a particular feature of the reciprocal width of the 
wavepacket, not shared by other parameters as the area 
of the soliton or its energy. We may find a physical ex- 
planation of this effect in equations iQ and ©: from 
the previous discussion is natural to assume that the in- 
tensity of radiation depends on the relative importance 
of the higher order dispersion terms (third derivative) 
to the internal kinetic energy of the wavepacket (second 
derivative). Since the ratio between the third order and 
second order dispersion terms is \2m*a3Ak/h'^\, a nar- 
rower wavefunction in the reciprocal space will radiate 
at a slower rate than a wider one. The limiting width, 
therefore, can be considered as a situation in which ra- 
diation takes place at a asymptotic slow rate. We have 
found evidence of this in our calculations, which show 
that the final wavefunction has a minimum value of the 
dispersion ratio. FigureOshows these ratios as a function 



of time, for the same parameters as in Fig. (21 The ratios 
have been calculated from the band structure according 
to the value of the effective mass and at the instanta- 
neous average momentum (shown in Fig^J and according 
to the instantaneous value of the reciprocal wavepacket 
width (Fig. EJ. 

As pointed out above, one of the consequences of the 
presence of higher order dispersion terms is to prevent the 
integrability of the system. As a result, the self-trapped 
solutions do not correspond to solitons but to the less 
restrictive class of solitary waves. One of the strict test 
of the soliton nature is found in the elasticity of their col- 
lisions. As a consequence, solitons are known to preserve 
their shape under such events. In order to explore this 
dynamical aspect, we have taken advantage of the ring 
scheme to model a collision experiment. In such geome- 
try, the transference of distinct momenta to wavepackets 
located at different regions of the waveguide seems fea- 
sible, simply by tailoring the proper phase mask for the 
phase imprinting. Also, the periodic boundary allows 
the wavepacket to decay into a solitary wave in a small 
spatial region, which is helpful when computing numer- 
ically problems in dimensions higher than one. In Fig. 
(71 we show the time evolution of two colliding solitary 
waves. The initial condition of the computation is a con- 
densate splitted into two wavepackets, at symmetric sites 
in the ring, and imprinted with opposite momenta. Each 
of the wavepackets resembles in number of atoms and 
initial momentum with the case shown in Fig. (21 Follow- 
ing the previous results of this paper, both wavepackets 
will decay into solitary waves at times close to 120ms. 
Therefore, the first plot in Fig. \7\ corresponds already 
to the two solitary waves in the ring heading one into 
each other. As it is apparent from the figure, the col- 
lision is highly inelastic, leading to the disintegration of 
the solitary wavepackets. 



IV. CONCLUSION 

We have analyzed the process of formation of solitary 
wavepackets of Bose-Einstein condensates loaded into op- 
tical lattices. Our results support the picture of matter 
radiation acting as a dissipative process that leads to the 
decay of the initial wavefunction into the solitary shape. 
We demonstrate that the source of radiation is linked to 
the presence of higher order dispersion terms inherent to 
the lattice band structure. Our point of view is consis- 
tent with previous studies on the effect of these dispersive 
terms in the propagation of light in optical fibers. After 
the process has been identified, we have demonstrated 
that the center of mass dynamics of the wavepacket is 
closely connected with the form of radiation through the 
conservation of the overall momentum, resulting in the 
acceleration or deceleration of the wavepacket as the soli- 
tary wave is formed. By radiating part of its mass, the 
wavepacket shrinks in reciprocal space, reducing accord- 
ingly the role of higher order dispersion terms. Conse- 
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quently, radiation diminishes in time as the wavepacket 
evolves to a sohtary wave of a fixed width in reciprocal 
space. This final width seems to be linked to the depth of 
the lattice sites and independent on the wavepacket pa- 
rameters (initial momentum and number of atoms). Fi- 
nally, we have studied the dynamics of collisions between 
these waves and found it highly inelastic, leading to the 
rapid dispersion of the solitary waves. All these results 
are supported by two-dimensional numerical integrations 
of the Gross-Pitaevskii equation in a ring-shaped optical 
lattice, and interpreted with the aid of a Wannier wave- 



function analysis. The size of the ring is chosen large 
enough to avoid particular effects related with this ge- 
ometry. 
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FIG. 1: Scheme for constructing ring-shaped optical lat- 
tices as proposed in the text: a laser (black arrow) is aimed 
perpendicularly to a ring-shaped atomic waveguide, after be- 
ing transmitted through an angular modulated transmission 
plate. In our calculations, the ring radius is 40/im, enclosing 
128 lattice sites. 



FIG. 2: (a) Plot of \ f{£,t)\'^ (wavefunction's envelope in the 
Wannier basis of the lowest energy band) as a function of 
time and the position in the ring, showing how the initial 
wavepacket converges into a solitary wave. The parts included 
in the dashed box have been magnified to enhance visibility, 
and the arrows show the radiation component traveling with 
the pulse, (b) Plot of the atom density at t = 180ms. The 
condensate of 100 atoms is initially released with momentum 
h58/ro along the ring circle, and the optical lattice is linearly 
built between t — 0.5ms and t — 10.5ms, afterwards remains 
constant with Vo/h ~ 128Hz. 



FIG. 3: (a) Plot of the evolution of the atom population in the 
first band of the Brillouin zone (projection into Bloch states 
as given by Eq. Q), for the same case as in Fig. |21 (b) The 
same but using only the part of the wavefunction ^ in the 
spatial box including only the region occupied by the solitary 
peak at every instant of time. Regions enclosed in the dashed 
box have been magnified. The reciprocal space momenta is 
given in units oi g = 2n/ao = 128/ro. 
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FIG. 4: Time evolution of the mean value of the reciprocal 
lattice momentum for the same parameters as in Fig. |21but 
initial momenta hko — h58/ro (solid line) and hko — TiQA/tq 
(i.e. at the edge of the Brillouin zone, dashed line). 



FIG. 5: Time evolution of the width of the wavefunction in 
the reciprocal space for different lattice depths: (a) Vb/fi — 
128ffz, (b) Vo/h ~ 2bQHz, (c)) Vo/h ~ UHz, and for con- 
densates with different initial conditions (number of atoms 
and momentum), written as {N,k/g} (a.l) {200,0.453}, 
(a.2) {100,0.453}, (a.3) {50,0.453}, (a.4) {100,0.5}, (b.l) 
{100,0.453}, (b.2) {50,0.453}, (c.l) {100,0.453}, (c.2) 
{50,0.453}. 



FIG. 6: Ratio between the third and second order dispersion 
for the same cases as in figure The ratios are computed 
at very instant of time using the widths shown in that fig- 
ure, and the corresponding effective mass computed from the 
instantaneous wavepacket's mean momentum. 



FIG. 7: Time evolution of two wavepackets with the same 
conditions as in Fig. |21 but opposite momenta, —h58/rQ and 
R58/ro. 
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